Reduction of optimum light power with Heisenberg-limited photon-counting noise in 

interferometric gravitational-wave detectors 
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We study how the behavior of quantum noise, presenting the fundamental limit on the sensitivity 
of interferometric gravitational- wave detectors, depends on properties of input states of light. We 
analyze the situation with specially prepared nonclassical input states which reduce the photon- 
counting noise to the Heisenberg limit. This results in a great reduction of the optimum light power 
needed to achieve the standard quantum limit, compared to the usual configuration. 
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Since the pioneering work by Caves it is well 

understood that two sources of quantum noise — the 
photon-counting noise and the radiation-pressure noise — 
constitute the fundamental limitation on the sensitiv- 
ity of an interferometric gravitational-wave detector. 
These limitations will be of potential importance in long- 
baseline interferometric detectors which are currently un- 
der construction (the LIGO project in the United 
States and the French-Italian VIRGO project |^ in Eu- 
rope are the largest ones). For example, the photon- 
counting shot noise will dominate at the gravitational- 
wave frequencies above 1 kHz in the VIRGO detector 
and above 200 Hz in the initial LIGO detector. With 
a further reduction of the thermal noise, planned in the 
advanced LIGO interferometer, the role of the shot noise 
will be even more important. 

For a coherent laser beam of light power P, the shot 
noise associated with photon-counting statistics scales as 
P~^/^ and the radiation-pressure noise scales as P^/^. 
The contributions of these two sources of noise will be 
equal for some optimum value Popt of the light power. 
Provided that classical sources of noise (such as thermal 
and seismic) are sufficiently suppressed, the interferom- 
eter with the optimum light power will work at the so- 
called standard quantum limit (SQL). A simple quantum 
calculation, based on the use of the Heisenberg uncer- 
tainty principle, gives the SQL for the measurement of 
the relative shift z = Z2 — zi in the positions of two end 



(Az)sQL = yj2hTlr 



(1) 



where m is the mass of each end mirror and t is the 
measurement time. For modern long-baseline interfer- 
ometers (like LIGO and VIRGO), Fabry-Perot cavities 
are used in the arms; so r is actually the cavity storage 
time, T = LJ- /nc, where L is the cavity length, J- is the 
finesse, and c is the velocity of light. The optimum power. 



for which the SQL is achieved, is 



P 



opt 



(2) 



where lo is the light angular frequency. For the initial 
LIGO configuration, the mirror mass is m ~ 11 kg, the 
cavity length is i ~ 4 km, the finesse is ^ ~ 200, and 
the wavelength of the Nd:YAG laser is A ~ 1.064 /im 
{uj ~ 1.77 X 10^*^ Hz). This gives the cavity storage time 
r ~ 8.5 X 10^^ s and an effective number of bounces 
b = TcjlL ~ 32. The corresponding optimum laser power 
is Popt — 191 kW and the SQL of the position shift mea- 
surement is (Az)sQL — 1-24 X 10~^^m. Achieving this 
SQL will make possible to measure gravitational waves 
with amplitudes h greater than ^ 3 x 10^^^. 

Presently, the available laser light power is insufficient 
for achieving the SQL (for example, in the initial LIGO 
configuration the input laser power is 6 W and the power 
recycling gain is about 30). Therefore, in advanced LIGO 
configurations, it is planned Q to reduce the shot noise 
by using more powerful lasers, in conjunction with the 
power-recycling technique 1^,0. However, for very high 
laser power, one encounters serious technical problems re- 
lated to nonuniform heating of the cavity mirrors caused 
by absorption of even a small portion of circulating light. 
The resulting thermal aberrations can seriously deterio- 
rate the performance of the interferometer |Q. Therefore, 
it will be very interesting to study possibilities for achiev- 
ing the SQL with low light power. 

The gravitational-wave detection community is quite 
familiar with the intriguing idea by Caves |^ to re- 
duce the photon-counting noise by squeezing the vacuum 
fluctuations at the unused input port. During the last 
decade, other interesting ideas has been developed in the 
field of theoretical quantum optics, based on the use of 
nonclassical photon states for the quantum noise reduc- 
tion in idealized optical interferometers |p|-p^ . The main 
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theoretical motivation of all those papers was to show the 
possibility of beating the shot-noise limit and achieving 
the fundamental Heisenberg limit for the photon-counting 
noise in an ideal interferometric measurement. 

The aim of the present work is to show that the op- 
timum light power needed for the SQL operation of an 
interferometric gravitational-wave detector with movable 
mirrors can be greatly reduced by the use of nonclassi- 
cal states of light with the Heisenberg-limited photon- 
counting noise. This result means that Heisenberg- 
limited interferometry is not only interesting for a demon- 
stration of the fundamental uncertainty principle, but can 
be also important for the experimental detection of grav- 
itational waves. 

Let us consider a long-baseline Michclson interferom- 
eter whose arms are equipped with high-finesse Fabry- 
Perot cavities, with end mirrors serving as free test 
masses. In the quantum description, two modes of the 
light field enter the interferometer through the two in- 
put ports of a 50-50 beam splitter. After being mixed 
in the beam splitter, the light modes spent time r in 
the Fabry-Perot cavities, and then leave the interferom- 
eter (through the same beam splitter, but in the oppo- 
site direction). The photons leaving the interferometer 
in the output modes are counted by two photodetectors. 
A gravitational wave incident on the interferometer will 
cause a relative shift z — Z2 — zi in the positions of two 
end mirrors, which results in the phase shift = (lut/L)z 
between the two arms. 

The performance of such an interferometer can be an- 
alyzed in the Heisenberg picture, using a nice group- 
theoretic description proposed by Yurke et al. Using 
the boson annihilation operators ai and a2 of the two 
input modes, one constructs the operators 



Jx = {a\a2 + 4ai)/2, 
Jy = -i{a\a2 - aJai)/2, 
Jz = {a\ai - ala2)/2. 



(3) 



These operators form the two-boson realization of the 
su(2) Lie algebra, [Jk,Ji] = i^UmJm- The Casimir oper- 
ator is a constant, — j{j + 1), for any unitary irre- 
ducible representation of the SU(2) group; so the repre- 
sentations are labeled by a single index j that takes the 
values j = 0, 1/2, 1, 3/2, . . .. The representation Hilbert 
space 7i| is spanned by the complete orthonormal basis 
I j, to) (to = j, j - 1, . . . , -j). Using Eq. (|), one finds 



^N{^N + 1), N = a\ai + 



ala2, 



(4) 



where N is the total number of photons entering the in- 
terferometer. We see that N is an SU(2) invariant; if the 
input state of the two- mode light field belongs to 7i| , then 
N = 2j. 

The actions of the interferometer elements on the 
column-vector J — {Jx, Jy, Jz)'^ can be represented as 



rotations in the 3-dimensional space [|. The first mixing 
in the beam splitter produces a rotation around the y axis 
by — 7r/2, with the transformation matrix Rj,(— 7r/2). The 
second mixing corresponds to the opposite rotation, with 
the transformation matrix Rj^(7r/2). The relative phase 
shift produces a rotation around the z axis by (f>, with the 
transformation matrix Rz {(p) . The overall transformation 
performed on J is the rotation by (j) around the x axis. 



R,(^) = R,(7r/2)R,(^)R,(-7r/2). 



(5) 



The information on the phase shift </> is inferred from 
the photon statistics of the output beams. Usually, one 
measures the difference between the number of photons 
in the two output modes. 



Qout — 



2 [(sin (t))Jy 



)Jz 



(6) 



If we assume that there are no losses in the interferome- 
ter and the classical sources of noise are well suppressed, 
then the uncertainty in the relative position shift z of the 
end mirrors is due to two factors The first one is 

the photon-counting noise. Indeed, since there are quan- 
tum fluctuations in gout, a phase shift is detectable only 
if it induces a change in (gout) which is larger than the 
uncertainty Agout- Consequently, the uncertainty in the 
phase shift due to the photon-counting noise is 



/ A ,\2 ^ (Agout)^ 
^ (9(gout)/a0)2- 



(7) 



If the detection is made on a dark fringe ((^ = 7r/2 in the 
unperturbed state), then the contribution of the photon- 
counting noise is 



(Az)^e 



A 



{Ajy 



pc 



{Jz)' 



Ape — I 



(8) 



The second source of noise is due to quantum fluctuations 
in the radiation pressure. The difference between the mo- 
menta transferred by light to the end mirrors, V = P2—Pi, 
is easily found to be 7^ = {2fiU)T / L)Jj.. The relative shift 
in the positions of the end mirrors, due to the transferred 
momenta, is (T/m)V. Therefore, the contribution of the 
radiation-pressure noise is 



(Az)2 = A,p(2AJ,)2, 



mL 



(9) 



Consider the standard case when the coherent laser 
beam of amplitude a enters the interferometer's one in- 
put port, while the vacuum enters the other. This in- 
put state I in) — |q!)i|0)2 (where |0) is the vacuum and 
|a) = exp(Q;a^ — Q;*a)|0) is the coherent state), satisfies 

(Jx) = {Jy) = O, (J,2) = (J2) = |a|V4, 
{Jz)^\a\y2, {N)=N^\a\\ 
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Using these resuhs, one finds 

(Az)2 = {Az)l, + {Az)% = Ap,N-' + ApN. (10) 
Optimizing (Az)^ as a function of N, one obtains 



(Az)2 ~ {ApcN-^ + A,pN) cosh2r. 



(14) 



N, 



(11) 



and Popt — fi^Nopt/T is given by Eg. (|2|), while the opti- 
mum value of Az is the SQL of Eq. dTfiT 

The characteristic noise behavior of Eq. (|l^) is some- 
times explained by the Poissonian photon statistics of the 
coherent state (i.e., by the fact that AA^i = {Niy^^, with 
Ni = a\ai). However, it is not difficult to see that this 
explanation is principally wrong. A simple calculation 
shows that if instead of the coherent state |a) at the first 
input port we will use an arbitrary state (pure or mix- 
ture) of the single-mode light field, the same result (0) 
will hold. This will be true, regardless of the statistical 
properties of the photon state at the first input port, as 
long as the vacuum enters the second input port. 

Caves proposed to reduce the optimum power 
needed to achieve the SQL by using the squeezed vac- 
uum in the second input port. If the carrier mode en- 
tering the first input port is in the coherent state \a), 
the two- mode input state is given by |in) = |a)i|f)2, 
where |0 = exp(ifat2 - ^^0^)10) (with ^ = re'") is 
the squeezed vacuum state. If one takes 6 = and real 
a, the input state satisfies 

(J.) = (Jy) = 0, (Jly) = [a^c^^^ + sinh^ r)/4, 
(Jz) = [a^ - sinh^ r)/2, N = + sinh^ r. 

In the usual situation, ^ sinh^ r, so one derives 



(Az) 



^pce 



-2r 



(12) 



The use of the squeezed vacuum reduces the photon- 
counting noise at the expense of the radiation-pressure 
noise. This results in the reduced optimum light power: 



Popt(r) ^Popt(r = 0)e 



-2r 



(13) 



while the optimum value of Az remains the SQL of 
Eq. (|l|). However, it is erroneous to think that the re- 
duction of the optimum light power is the merit of the 
squeezed vacuum alone; actually, the state of the car- 
rier mode is important as well. While the state in the 
carrier mode was mixed with the phase-insensitive vac- 
uum, the behavior of the quantum noise was determined 
by the mean number of carrier photons only. However, 
when mixing the carrier mode with a phase-sensitive state 
(e.g., with the squeezed vacuum), the precise matching 
between the quantum states of the two modes is impor- 
tant. For example, if the state of the carrier mode satisfies 
{a\^ + al) = and (A^i) sinh^ r, then we obtain 



This results in the same optimum power as for the 
normal vacuum case, but the sensitivity deteriorates: 
{Az)lp^ (Az)gQLC0sh2r. For example, this will be 
the case for the carrier mode in the coherent state \a) 
with arga = 7r/4 or in any phase-insensitive state (a 
state is called phase-insensitive, if its density matrix is 
diagonal in the Fock basis; examples are the Fock states 
themselves or the thermal state). On the other hand, 
if the carrier mode in an arbitrary state is mixed with a 
phase-insensitive state (or, more generally, with any state 
satisfying {a^) = 0), then 

(Az)2 = (2iViiV2 + Ni+ N2) [Ap,{Ni - N2)-^ + Ap] . 

Here, we used notation Nk = {a\ak), k = 1,2. Clearly, 
the quantum noise cannot be reduced here, compared to 
the vacuum case; in particular, for A'^i ^ N2, the op- 
timum power remains as in Eq. (^), but the sensitivity 
deteriorates: (Az)^pt ~ (A2)|ql(1 -I- 2N2). 

From the above arguments, one understands that the 
reduction of the optimum power can be achieved with 
a proper phase matching between the two input modes. 
In this relation, it is interesting to consider input states 
which lead to the Heisenberg-limited photon-counting 
noise [||-^3|_. It is well known that the shot-noise limit 
{A(/))pc = N"-^/^, achieved with the vacuum at the sec- 
ond input port, is not a fundamental one. Using the 
uncertainty relation {AJj;){AJy) > ^\{Jz)\, one obtains 
(A0)pc > (2AJa;)~-^. Since for any input state |in) S Hj 
the relation (AJ^)^ < + 1) holds, one finds the 

Heisenberg limit 



(A0)pc > [2j{j + 1)] 



-1/2 



(15) 



Consequently, for large photon numbers {N = 2j ^ 1), 
the photon-counting noise (Az)pc scales as 1/P. 

It can be shown that the shot-noise limit can 
be surpassed with the so-called intelligent (minimum- 
uncertainty) states ||l^,|l^. (The use of intelligent 
states for achieving the Heisenberg limit in spectroscopy 
was discussed in [Q.) The Jx-Jy intelligent states, 
by their definition, equalize the uncertainty relation: 
{AJx){AJy) = ||(>/z)|- These states are determined by 
the eigenvalue equation 



{ijJx - iJy)\X,r]) = A|A,77). 



(16) 



The spectrum is discrete: A = imo-\/l — 77^, where 
'Ti'O = Jt] ~ li ■ ■ • J a-iid 77 is a real parameter given 
by I?/! = AJy/AJx- Recently, a method for experimental 
generation of the SU(2) intelligent states was proposed in 
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If the two-mode light field entering the interferometer 
is prepared in the Jx-Jy intelligent state, the quantum 
noise takes the form 



(Az)2 = Ape(2A J, 



Ap(2AJ,)2 



(17) 



For \ri\ < 1, the intelligent states are squeezed in Jy 
and anti-squeezed in J^, thereby reducing the photon- 
counting noise below the shot- noise limit, on account of 
increasing contribution of the radiation-pressure noise. 
For 77 ^ 0, one obtains [M 



(2AJ,)2 = 2|(J,)/77| ~ 2if -ml+ j). 



(18) 



and the Heisenberg limit for the photon-counting noise is 
achieved when toq = 0. Then, for large photon numbers 
(N = 2j > 1), we obtain 



(Az) 



2Ap,N~ 



(19) 



Optimizing (Az)^ as a function of N, we find (Az)opt — 
(Az)sQL, but the optimum light power needed to achieve 
the SQL is dramatically reduced: 



/2mL2 



1/2 



-fopt — 



2\ 1/2 



(20) 



Using the parameters for the initial LIGO, we find the val- 
ues iVopt ^ 4.3 X 10^" and Popt — 9 fiW. Compare these 
values with those obtained in the standard configuration 
(with the vacuum at the second port): iVopt — 9.2 x 10^" 
and Popt — 191 kW. We see that the use of the intelli- 
gent states with the Heisenberg-limited photon-counting 
noise can in principle reduce the optimum light power by 
a factor ~ 2 x 10^°. 

There were proposals to achieve the Heisenberg limit 
for the photon-counting noise by driving the interferom- 
eter with two Fock states containing equal numbers of 
photons |^,|l3|. (The use of Fock states was also pro- 
posed for Heisenberg-limited interferometry with matter 
waves 1 17 and for Heisenberg-limited spectroscopy with 
degenerate Bose-Einstein gases [|l^.) The corresponding 
input state is |in) = |n)i|n)2 = |i, 0) with j — n — N/2. 
Clearly, this input state cannot be used when one mea- 
sures the photon difference gout at the output. However, 
as was shown in Ref. [ p^ , the Heisenberg-limited photon- 
counting noise is achieved by measuring the squared dif- 
ference S — Qout = 4J^out the output. (We do not 
discuss here technical problems involved in this kind of 
measurement.) In this case, the uncertainty in the phase 
shift due to the photon-counting noise is 



tan 



2 



2 — tan 
4j(.? + 1) 



(21) 



For (j) — Q (this corresponds to a dark fringe for the 
measurement of S), the Heisenberg limit is achieved: 



(A0)pc = [2j(j -I- 1)]~^/^. Of course, this improvement is 
on account of the corresponding increase in the radiation- 
pressure noise, because (AJ^,)^ = \ + takes its max- 
imum value. Therefore, for large photon numbers, we ob- 
tain (Az)po - {2Ap,YI^/N and (Az),.p ~ (Arp/2)i/27V, 
recovering the result of Eq. (p^. In fact, this quantum 
noise behavior and the corresponding reduction of the op- 
timum power are characteristic for input states with the 
Heisenberg-limited photon-counting noise. 

It should be emphasized that above results are valid for 
a lossless interferometer, while mirrors of realistic detec- 
tors (for example, LIGO) do have losses. The Heisenberg- 
limited photon-counting noise can be achieved only if the 
losses are sufficiently small. Let F be the dimensionless 
coefficient of losses, defined by TVout = Ne~^ ~ N{l — T). 
In the case of input Fock states, a simple analysis shows 
Ip^ that the Heisenberg-limited photon-counting noise 
(A(?!))pc N-'^ can be obtained only for iVF < 1/2. 
Sure, the value F ~ 10~^^ is impossible to achieve with 
the present technology. The problem of losses is of great 
practical importance but nevertheless it does not cross 
out the principal value of the idea to reduce the opti- 
mum light power by using nonclassical input states. For 
example, one can imagine a realization of this idea with 
a prototype interferometer which should have smaller m 
and L and larger r. 

In conclusion, we analyzed the behavior of quan- 
tum noise, which limits the sensitivity of interferomet- 
ric gravitational- wave detectors, for various input states 
of the light field. We found that by using nonclassical 
input states exhibiting the Heisenberg-limited photon- 
counting noise, the optimum light power needed to 
achieve the SQL can be significantly reduced, compared 
to the usual configuration. Of course, a practical real- 
ization of Heisenberg-limited interferometry will depend 
on future theoretical and experimental progress in the 
methods for production of stable and sufficiently power- 
ful sources of nonclassical light and on technical solutions 
for the reduction of losses and classical sources of noise. 
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